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INTRODUCTION 

In the present article we investigate some geophysical aspects of the tidal 
perturbations in the motion of artificial satellites and develop a system of formulas 
convenient for computation of the tidal effects in the elements using a step-by- 
step numerical integration. 

Recently the tidal perturbations in the motion of artificial satellites captured 
the imagination of theoreticians and observers. The significance of these pertur- 
bations lies in the fact that they are intimately connected with the elastic properties 
of the Earth and with the distribution of density inside the Earth. Thus, the 
comparison between the observed and computed tidal perturbations of the satellite 
permits one to obtain information about the average elastic properties of the 
Earth and to check the assumed model of the Earth. 

The attraction of the Earth by the Moon and the Sun produces elastic dis- 
placements inside the Earth body {bodily tides) and the tidal buldge on its surface 
(surface tides). The result is a small change in the distribution of mass and, 
consequently, small variations in the exterior geopotential. These variations 
constitute the exterior perturbing tidal potential, and they produce the perturba- 
tions in the motion of an artificial satellite. 
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The tidal perturbations in the motion of the satellite are not the result of the 


gravitational attraction by the surface buldge only. The bodily tides contribute 
their share of influence {approximately 40%). 

The exterior tidal potential can be expanded into a series of products of 
spherical harmonics. The first factor in each product is a polynomial in the 
components of the unit vector directed from the center of the Earth toward the 
Moon (the Sun). The second factor, of the same degree and order as the first, 
is a similar pol 5 momial in the components of the unit vector directed toward 
the satellite. This form of the tidal potential is convenient for the computation 
of perturbations by means of numerical step by step integration. If we prefer 
an anal 3 d;ical expansion, then, to facilitate the integration, the products of 
spherical harmonics are expanded into Fourier series with arguments which 
are linear combinations of the arguments^ , I', F, D and r of the lunar theory 
and of the mean ai^lar elements of the satellite (Musen and Estes, 1971). All 
these arguments are either linear or nearly linear in time and, consequently, 
the integration in analytical form does not present any difficulty. The program 
for manipulating the Fourier series which appear in the theory of tidal perturba- 
tions was developed by C. Hipkins and R. Estes of BTS Company. 

Each term of the Fourier expansion contains a factor, i.e. the Love number, 
which represents a measure of the elastic response of the Earth to the given 
tidal frequency. These numbers depend on the laws of change of Lame elastic 
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parameters and the density inside the Earth, on the degree of the spherical 
harmonics, and,to a lesser extent, on its order and on the frequency itself 
{Alterman et al.» 1959). 

EXPANSION- OF THE EXTERIOR TIDAL POTENTIAL 

The tidal oscillations of the point r " in the earth are governed by the partial 
differential equation (Alterman et al., 1959): 


32 s 

p ;; pl^ + gV- (ps) - V(ps * g) + V’ a, 

3t2 


( 1 ) 


where V is the del-operator relative to r", a- is the stress tensor 

err: ja(Vs + sV) + MV-s, 

I is the idemfactor, g is the undisturbed acceleration of gravity, and is the 
sum of the direct tidal potential acting on r" and the potential due to the tidal 
disturbance in the interior geopotential, as caused by the presence of elastic 
displacements inside the Earth and on its surface. The density p and the Lame 
elastic parameters, \ and m, are considered in the frame of the present work as 
functions of v" - r"/R only, where R is the mean radius of the Earth. The po- 
tential satisfies Poisson’s equation 

^ + 4t 7GV- (ps). (3) 

The particular solution of (3) 


n 


- G f dv. 

Ju, I r - r " 1 


(4) 
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represents the exterior tidal potential acting on the satellite, where r is the 
position vector of the satellite. The integral in (4) is taken over the volume 
of the whole tidally disturbed Earth. In our approximation we assume that the 
tidally undisturbed Earth is a sphere of the radius R and neglect the influence of 
the Coriolis force. This assumption, however, must be removed in the investi- 
gations about the response of the satellite to the tidal forces over a very long 
interval of time. Under these assumptions the potential (4) represents the effects 
of the body and surface "solid" Earth tides on the motion of the satellite. The 
effects of the oceanic tides are not considered in the present paper. 

We split the integral (4) into the sum of the integral over the volume of the 
sphere r" = R and the integral over the surfaces of density discontinuities, 


n = - G 


V*(yOS) , ^ 

1 dv + G 


r - r 


f p's- r"0 

Js It- r"| 


dS. 


( 5 ) 


where r"° is the unit vector in the direction of r", 

dv = r "2 dr "dr", 
dS =r"^dcr", 
do-" = cos ,e"dyS"d^", 

and/3 is the latitude and A." the east longitude^ and p’ is the density discontinuity, 
p' = P on the outer surface. 

The first term in (5) represents the disturbing influence of the internal 
elastic tidal displacements inside the Earth's body; the second term represents 
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the perturbative influence of the tidal buldge. The displacement s , as given by 
the differential equation (1), is a combination of the spheroidal and torsional 
oscillations of the Earth. 

In the investigations of the influence of tides on the motion of a satellite it 
is convenient to represent the components of the spheroidal oscillations in the 
form: 

s (K) = _i!!_ — ^(v") Y„„(yS", K") r"° 


( 6 ) 


r"L„,(v") gradY„J/3", \“)} k^Jk) 


This representation of differs in form only slightly from the corresponding 
representation by Alterman et al (1959). In (6) m'/M is the ratio of the mass of 
the Moon (Sun) to the mass of the Earth, and p' is the lunar (solar) parallactic 
factor R/a', where a' is the mean distance of the Moon (Sun) from the Earth. 

H (i'") and L (v") can be termed the generalized Love and Shida "numbers", 
with u" = r"/R. They satisfy a system of ordinary differential equations (Alterman, 
et al., 1959), (Takeuchi, 1950) which can be integrated only numerically. The 
functions Y^^ (/9", \") are the normalized spherical harmonics, 






2 n + 1 , (n - |m| ) ! 

(n -f- I m| ) ! 


Pjj(sin yS") e 


„ + i m\ . 


(7) 
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p;(x) (1 - - 1)" 

n ! 2" dx"'^'" 


with the normalization and orthogonality conditions 


r2n r+TT/2 

J “""J 

•'O •'-77/2 


.(/S". ^“) cos /3"d/3" = 8 ,8 


( 8 ) 


In spherical coordinates we have: 


= + - 
nm 


3Y, 


9Y 


( 9 ) 


o/?“ r“cos/3" 


where e^„ and e^„ are the corresponding unit vectors. The factor . 
exp i( /ft+m^in (6) is inherited by in {6) from the "lunar" ("solar") factors (a'/r')’^ 
^nra in the expansion of the static tidal potential acting on each point of 

the Earth, The part is a linear combination of -i', F, D and r, and 6 is the 
sidereal time. Besides (6) there is a second geophysically important solution 
of (1). It is of the form 


s = 

nm 


4tt 

2n + 1 


Rp'n+I^ (r" X grady 


and represents the torsional oscillations of the Earth (Alterman et al., 1959), 
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Jeffreys (1967). However, it is not of any importance in ou'r problem. Its radial 
component is zero and, because ’ ’ ■ 

div(W„,r"x gradY„J = 0. 

it does not produce any, dilatation. Consequently, it will not produce any dis- 
turbances in the exterior geopotential at all, and will not influence the motion 
of an artificial satellite. Thus, only the spheroidal osciUatibhs of'the iJarth as ■ 
given by (6) can produce perturbations in the motion of' a^satellite. ri*- ' 

The tidal potential (5) can be represented as a sum of particular tidal po- 
tentials of the form: 


n = - G 

nm 


|r - r"| 


dv + G 


-JHH dS. 

Js |r-r"I 


( 10 ) 


We have 


V-(ps ) = - ’+ pV«, 

R dy" 


(U) 


From the last equation, taking into account 


r''0.VY„„ = 0> 


. ,<!• 


V^Y = - 


+ 1) Y 

«2 


and (6), we obtain: 
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V- s = — • iL p'"+i 
2n + 1 M 


I 9 

— 


_ T V //?" \ "^ i»i(^t+iii0) 

« “*^0111 nm ’ * 


Substituting (12) into (11), expanding | r” - r | • in terms of spherical harmonics, 
and taking (6) into consideration we obtain: 

n = (pp'r^ ei'<‘ K ra 

^ 2n + 1 R “ 


. ^ (^) 
• nm ^ 


Y_ (S, a), 




where k is the Love number. The summation symbol relates to the surfaces 

nm 

of the density discontinuity and we set 




L„.(V) t 

/I nm ' ' nm . « 

v" dv 
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p is the parallactic factor R/a, and a is the semi-major axis of the satellite’s 
orbit, is the mean density of the Earth. The numerical value of the first 
term in (14) gives a general idea about the order of magnitude of the perturbing 
by the tidal buldge alone. The value of the second term gives information about 
the magnitude of the disturbing influence of bodily tides. In (13) all short period 
terms depending on the Earth’s rotation and involving the sidereal time have 
disappeared automatically. 

The most significant tidal perturbations in the motion of a satellite are those 
of long period. We obtain the disturbing potentials with long period terms 

only, by averaging the expressions (13) over the orbit of the satellite. Thus: 


2n + 1 R ^ ^ 


k e 

nm 


iKt 


w 


where 


W 

nm 



a) dg 


(17) 


and g is the mean anomaly of the satellite. The results of this averaging in 
terms of the orbital elements of the satellite are given by several authors 
(Kaula, 1969), (Kozai, 1965), (Newton, 1968), (Musen and Estes, 1972), (Musen 
and Felsentreger, 1973) and therefore are omitted in the present exposition. 
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The differential equation (1) indicates the dependence of Love numbers 
on the tidal frequency k. 

The dependence of on n is much stronger than the dependence on m and 
K. As a consequence the satellite will have some "difficulties’* distinguishing 
between different k^^ for a given n. For this reason static model of the Earth 
is being used at the present time in the computation of the tidal effects in the 
motion of satellites and the dependence of Love numbers on the tidal frequencies 
is suppressed. 

At present the value k is used for all m and all k associated with the 
Fourier expansions of 



a') 


/a'X" 

md _) Y* (S', a'). 
\r V 


in terms of the arguments , 't' , F, D and r of the lunar theory. Under these 
limitations, at the present time the disturbing function of the form 


Q 


4tt 

2n + 1 


R 


(13*) 








is being used, 
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and we deduce for the total tidal potential 


n 


+ 00 111= + n 

■zz: 

n =2 




( 18 ) 


by taking the addition theorem 


m^ + n 


S) = 2^ y;.(S-, .') y„.(5. a). 


(19) 


cos S = r*’ ■ T ' ® , 


into consideration. 


m= 2 


( 20 ) 


The expansion (20) is the standard one and presently serves as the basis for the 

1 

computation of the tidal perturbations in the motion of the artificial satellite as 
caused by the tides of the "solid Earth". 

However, with improved observational techniques and the use of satellite 
altimetry tiie dependence of Love numbers on frequency must be considered 
for accurate geoid determination from satellite observations. In the present 
exposition we develop a set of differential equations for the long period tidal 
perturbations in satellite elements in a form suitable for numerical step-by-step 
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integration. Having this goal in mind we re-write the disturbing function (13) as 


n 3 

2n + 1 




(r 


Y„„(S, a). 


( 21 ) 


where we set 

O (/<)= — Ik K 

and take the relation 

GM = j»2a3 

into account. 

If for given n the Love numbers can be replaced by their average values , 
then instead of (21) we can use the "combined" disturbing function (13’), which 
we re-write in the form: 


2n + 1 


n^aA'* 

nm 


Y„.(S. a). 


(22) 


where 




q =— Ff’p'"+lk . 
^ M 


( 23 ) 
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Representations (21) and (22) lead to a conapact form of the differential equations 
for the perturbations in elements. 

We assume that the spherical harmonics y (S, a), y* (S', a') are repre- 

nm nm 

sented as polynomials in the equatorial components 

K - cos 8 cos a, fji = cos 8 sin a, = sin S, 

= cos S' cos a', /x' = cos S' sin a', v' sin S' 

of the unit-vectors r° and r'® , respectively. 

DIFFERENTIAL EQUATIONS FOR THE VARIATION OF ELEMENTS 

At the present time only the long period tidal perturbations in elements 
can be easily observed. It is convenient to compute them by means of a 
step-by-step numerical integration making use of the Gaussian form of the 
differential equations for the variation of elliptic elements: 

de _ /l - . 1 (24) 

d t fi a e 277) a 


drF 

dt 



S cos f + 



T s in 


dg 


(25) 


+ 2 s i n^ 


1 

2 


da 

dt 
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(26) 


dt 


1 

-}- ■ — 

7ia/l - sin i 


1 C’ 

2" Jo = 


2 sin(f + cl>) dg, 


di 

dT 


1 

H — • — 

n a/l - 


2tt 



- Z cos(f + co) dg 
a 


(27) 


dL 

dt 


n a 



(28) 


+ 2/1 - e2 


s in^ 


1 dfl 

2 dt ’ 


where 

S = r° F, T= n®F, Z = BF 
are the projections of the disturbing force 


(29) 


F = Vn 

on the directions 

r° = Pcos f + Qsin f, n° = Rx r® = -Psin f + Qcos f, and R, 

P , 0 and R are the Gibbsian vectors associated with the instantaneous set of 
elliptic elements. The values of the components S, T and Z are computed for 
a set of points along the instantaneous orbit as defined by this set. 
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We have: 


V = rO ^ ^ 


3r 


-(I-rOrO)*Vo, 


(31) 


where is the del-operator relative to r® , 



+ j 


3 , B 


Making use of (21) we obtain for the disturbing forces 


F = vn ^ 

ntn ncn -|- 1 




n+2 

{- (n + 1) r° + (I - r®r°)- Vg} (^^) 


We deduce from (24) — (28), (30) and (32) a set of formulas convenient for 
programming: 


de 

dt 


^TTn 

Jn + 1 


1 

e 







B 

nm 


df , 


(33) 


d?7 _ 477'ra 

d?" 2n + 1 


1 q • e ^ ^ ‘ 

A 


(34) 


|*A„,cosf + 



B„„ sin f ^ df + 2 s in^ _ 


d£l 
dt ’ 
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1 


d£l _ 477-n 

dt 2n + 1 


(1 - e^) sin i 




( 35 ) 


2n 


f 2rr 

^nn. + ^ df, 


d i _ 477'n 

dt “ 2n + 1 





iKt 


(36) 


. 1 

2tt 



COS(f + 77 


a) df. 


dt 


47771 

+ 

2n + 1 



"Inm • e 


iKt 


1 

2t7 j 


riir 


D 

nm 


df 


(37) 


+ 


e 

1 + /F- e2 



+ sin^l 

2 dt 
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where we set for brevity; 


^nn. = + 1) 



Y , 

nm 


and 


(fj ‘ 

D = Y 

I j- / nm 

' Vq - R X 

= (.-R^ - „RJ ^ ♦ (>.R_ - ^ 


.frR,-^R^)|_, 

If instead of particular values of Love parameters (/<) we decide to use the 
average value , which assumes k^^ («r) to be the same for a given n and for 
all admissible m and k , then we can also use a more compact form of the 
disturbing function: 
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( 22 ) 


n = 


4 tt 
2n + 1 


rt%A* * 

nm 


( r ‘ 


where 



Y' 

nm 


(23) 


% = 


M 


k 


n 


and the differential equations (33) - (37) become: 


de 

dt 


4tt 

2n + 1 


• — A' * — 
e 2 tt 



B df, 

nm ’ 


(33’) 


^ = • 1a'* 

dt 2n + 1 e 


1. r 
2-1 


f + L + ] sin t i-df 


1 - eV 


in t^c 


j ■ 2 d 

+ I sin-'— , 

2 dt 
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n 


(35') 


_ ^ 4t7 

" ^ 2n + 1 


(1 - e^) sin 


A'* 

nm 


27T 



C 


nm 


sin(f + 77 - A) df, 


^- 477 . « 1 ( 36 >) 

dt 2n + 1 , 2 * 277 I + 77^ - <n.) df, 

i - e Jq 


dL 

dt 


+ 


4t7 

2n + 1 


’ A ’ * ■ ^ 

/l _ g2 2 t7 



df 


(37') 


+ 


e 

1 + / I - e2 



+ 2/1 - e2 sin^i 

2 dt 


Conclusions 

In expanding the exterior tidal potential in the present work we assumed 
spherical symmetry of the Earth’s elastic properties and density. The Love 
number for a given frequency can be computed (Alterman et al., 1957) on the basis of 
the assumed laws (from seismic evidence) of variation of Lame elastic parameters, 
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A. and /u., and of the density with depth. At the moment these hypotheses, commonly 
used in the theory of oscillations of the Earth, seem to be plausible. 

At the present time, however, an even more simplified model of the exterior 
tidal potential, which is implicitly based on a static model of the Earth, is being 
widely used in the computations of satellite tidal perturbations. In the frame of 
this model the dependence of Love numbers on the order of the spherical har- 
monics and frequency is supressed in the expansion of the tidal potential and 
the same Love number is assigned to all spherical harmonics of the same degree. 

With the accumulation of a long series of observations and a further im- 
provement in observational techniques, consideration should be given to the 
dependence of Love numbers on frequency, at least for the few most important 
lunar and solar tidal constituents. 

Only the investigations based on the use of an analytical (or semi-analytical) 
theory can easily provide information on the significance of this dependence. 

These investigations of the dependence of Love numbers on frequencies are 
necessary in connection with the determination of an accurate geoid. 

There are other un-resolved problems which await solution. The influence 
of the core is one such problem. The speed of rotation of the core is different 
from the speed of rotation of the mantle and crust. 

As a consequence resonances and variability of amplitudes and phases 
appear in the expansion of the tidal effects in the motion of a satellite. The 
determination of lags of tidal constituents in the perturbations satellites will 
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permit one to xmderstand the process of dissipation of energy. The problem of the 
influences of oceanic tides on satellite motion still is not solved numerically, 
although it is satisfactorily understood theoretically. The numerical solution 
of this problem is tied with the integration of the Laplace tidal differential 
equations over the global ocean for different tidal constituents. 

Our final conclusion is that in parallel with attempts to determine the 
average elastic properties of the Earth from satellite observations it is neces- 
sary to approach the problem from the geophysical side and integrate the dif- 
ferential equations for Love parameters, at least for the most important tidal 
frequencies in the motion of the satellite. 

These are some problems which constitute topics for a long term future 
work. 
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